
  x = [0  .999999999999934  -.499999019114416  -.500002802182133  1.54685516241072e-006  ];

  y = [0  0.  -.866025970098176  .866023785937356  -1.39078549577776e-007  ];


  z = [1  3.64481743070517e-007  -3.68275833545306e-007  -8.62773940895901e-007  -.999999999998794  ];
  
   plot3(x,y,z,'rs','Markersize',15);
  sphere;
  axis equal;
  hold on;
  plot3(x,y,z,'rs','Markersize',15,'Markerfacecolor','black');
  hold off;
  
  s = 0;
  
  for i = 1:length(x)
      for j = i+1:length(x)
          s = s + (x(i)-x(j))^2 + (y(i)-y(j))^2 + (z(i)-z(j))^2;
      end
  end
  s